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Abstract. We investigate the possibility that radio halos 
of clusters of galaxies are caused by synchrotron emission 
of cosmic ray electrons (CRe), which were produced by 
cosmic ray protons (CRp) interacting hadronically with 
the intra-cluster medium (ICM) protons. We perform cos- 



vannini ( [L996|) Fer etti fll999| ), Enfilin (|l999b|) , and Gio- 
vannini et al. (|199S| ). 

Whereas the cluster relic sources are nowadays widely 
believed to trace the presence of large shock waves in the 
ICM flEnfilin et al., 1998] |Roettiger et al., 1999] ; |Enfilin & 



obtain a sample of ten magnetized galaxy clusters. They 
provide realistic models of the gas and magnetic field dis- 
tribution, needed to predict the CRe production rates, 
their cooling, and their synchrotron emissivity. We as- 
sume a CRp population within the ICM with an energy 
density which has a constant ratio to thermal energy den- 
sity. This ratio is adjusted in such a way that one of the 
simulated clusters reproduces the radio luminosity of the 
radio halo of the Coma cluster of galaxies. Our model ex- 
hibits the observed low degree of radio polarization and 



mologi[;al magneto- hydrodynamics (MHD) simulations to Krishna, 2000[) the origin of radio halos is less clear. The 



lifetime of the radio emitting CRe is much shorter than 
any reasonable transport time over cluster scales. This 
means that the CRe need to be continuously re-energized, 
or freshly injected on a cluster- wide scale. The latter would 
happen if the CRe result from the decay of charged pions 
produced in hadronic interaction of energetic CRp dis- 
tributed throughout the cluster volume ( Dcnnison, 1980; 
Vcstrand, 198 2|; |Valtaoja, 1984| ; |Schlickeiser et al., 1987 ; 
Blasi & Colafrancesco, 1999). The CRp have lifetimes in 



the ICM of the order of a Hubble time ( Berezinsky et al 



We provide estimates for the expected gamma ray and 
neutrino flux. The necessary CRp/thermal energy ratio 
is 4 ... 14 % (£p,,„i„/GeVr - 375 (for the range of mag- 
netic field strengths suggested by Faraday measurements) , 
where E^mm is the lower kinetic energy cutoff of the CRp 
with spectral index a p ~ 2.375. Assuming this ratio to be 
the same in the whole set of simulated clusters a T x — L u 
relation is predicted which follows the observed relation 
well. 



Key words: Magnetic Fields - MHD - Radiation mcch- 
anism: 



has a similar radial emission profile as the Coma cluster. 1997; Enfilin et al., 1997), giving them enough time to dif- 



fuse away from their source and to maintain a cluster wide 
distribution. 

It is the purpose of this work to demonstrate that 
a simple model for hadronic electron injection in a re- 
alistic magnetic field configuration leads to radio halos 
which reproduce several observations: the profile of the 
radialy decreasing radio emission, the low radio polariza- 
tion, the radio to X-ray surface brightness correlation, and 
the cluster radio halo luminosity-temperature relation. We 
do not attempt to reproduce all observational information 
on cluster radio halos. For example the observed spatial 
spectral index variations of individual clusters ( piovan- 



non-thcrmal - Galaxies: intergalactic medium nini et al., 1993| ; peiss et al., 1997[ ) are not reproduced by 



Galaxies: cluster: general - Radio continuum: general 



1. Introduction 

Some clusters of galaxies contain large, Mpc-sized regions 
of diffuse radio emission, the so called cluster radio halos 
and cluster radio relics. Although known for decades, these 
sources receive increasing attention due to the improved 
detection sensitivities and their important theoretical im- 
plications for the no n-ther mal con tent of the ICM. For 
reviews see Valtaoja (|l984| ), Jaffe ( |l992[ ), Feretti & Gio- 



the presented model. It is our goal to investigate the ob- 
servational consequences of a conceptually simple model. 
Sophisticated modifications, which fine-tune the model to 
the data, are principally possible. But first an understand- 
ing of idealized models is most helpful. The assumptions 
of our model are: 

1. CRe are secondary particles of hadronic interactions of 
CRp with the background protons of the ICM. The ori- 
gin of the CRp is not specified here, but a likely expla- 
nation is particle acceleration in shock waves of clus- 
ter merger events and in accretion shocks of the clus- 
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ters flColafrancesco fc Blasi, 1998| ). ICM CRp might 
also migrate from radio plasma of radio galaxies into 



the IC M flValtaoja, 1984|; [Enfilin et al., 1997| ; |Enfilin : 
199|9a|; Blasi fc Colafrancesco, 1999|) or result fro m su- 
pernova driven galactic winds (Volk et al., 1996). 



2. The energy density in CRp is assumed to be propor- 
tional to the energy density of the thermal ICM gas. 
This is reasonable if the thermal gas and the CRp were 
energized in the same shock waves of cluster mergers, 
since a constant fraction of energy dissipated in the 
shock should go into the CRp population. Previous 



work ignore d the spatial profiles ( Dennison, 1980 ; Vcs- 



trand, 1982 ) or assumed one single central point source 



from which the CRp diffused away, resulting in a cen- 



trally very peaked CRp profile flBlasi 
1999| ). 



Colafrancesco, 



meon 
min/max 
Sim. 9 



The CRp spectral index is independent of position or 
cluster and therefore the radio spectral index is con- 
stant in this model. Although this is certainly an over- 
simplification since e.g. the radio halo of the Coma 
cluster does show spectral index variations, it allows 
to predict radio halo statistics. 

No significant re-acceleration or particle diffusion of 
the CRe is occurring. 

The radio spectrum is assumed to be quasi-stationary, 
since the electron cooling time is short and a stationary 
CRe population is rapidly established. The spectrum 
of this population results from the shape of the pro- 
ton spectrum, which is assumed to be a single power- 
law, and the synchrotron and inverse Compton cooling. 
Hence electron and radio spectra are also power laws. 
The magnetic field configuration results from amplifi- 
cation of pr e-cluster seed fields as simulated by Dolag 
et al. ( 1999 ). We therefore investigate very inhomoge- 
neously magnetised clusters. Although the importance 
of realistic field configurations for the a ppearence of 
radio halos was stressed by Tribble ( 1991 ), even recent 
work on halos of hadronic origin assumed spatially con- 
stant field strength. 



2. MHD Cluster Formation Simulation 

We used the cosmological MHD code described in Dolag et 
al. (1999) to simulate the formation of magnetised galaxy 
clusters from an initial density perturbation field. The 
evolution of the magnetic field is followed from an ini- 
tial seed field. This field is amplified by the compression 
during the cluster collapse. Merger events and sheer flows, 
which are very common in the cosmological environment 
of large-scale structure, lead to Kelvin-Hclmholtz instabil- 
ities. They further increase the field strength by a large 
factor. In order to reproduce the present fiG fields an ini- 
tial nG field strength is required. 
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Fig. 1. The radial profile of the magnetic field strength 
of the simulated clusters is shown in units of the virial 
radius i?200- The solid line represents the field strength 
in the cluster Sim. 9. The thick dotted line is the profile 
averaged over all ten simulations, the thin dotted lines are 
the minimum and maximum values found in the different 
simulations. 



2.1. GrapeSPH 

The code combines the merely gravitational interaction of 
a dark-matter component with the hydrodynamics of a 
gaseous component. The gravitational interaction of the 
particles is evaluated on GRAPE boards, while the gas 
dynamics is computed in the SPH approximation. It was 
also supplemented with the magneto-hydrodynamic equa- 
tions to trace the evolution of the magnetic fields which 
are frozen into the motion of the gas because of its as- 
sumed ideal electric conductivity. The back-reaction of the 
magnetic field on the gas is included. Extensive tests of the 
code were successfully performed and described in the pre- 
vious paper. V • B is always negligible compared to the 
magnetic field divided by a typical length scale. The code 
also assumes the ICM to be an ideal gas with an adiabatic 
index of 7 = 5/3 and neglects cooling. The surroundings 
of the clusters are dynamically important because of tidal 
influences and the details of the merger history. In order 
to account for this the cluster simulation volumes are sur- 
rounded by a layer of boundary SPH particles in order to 
represent accurately the sources of the tidal fields in the 
cluster neighborhood. For details of the code, the models 
and the obtained magnetic field structure see Dolag et al. 
(1999; in preparation). 
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Fig. 2. Comparison of synthetic Faraday-rotation mea- 
surements with a sample of measurements in Abell clus- 
ters. The absolute values of Faraday rotation measure- 



ments (Kim et al., 1991) vs. radius in units of the Abell 
radius are shown. Obviously, the dispersion increases to- 
wards the cluster center. The solid curves mark the me- 
dian of the measurements, the dotted curve is the median 
obtained from the simulated cluster sample. Some mea- 
surements have Faraday measurements which are outside 
the region shown (marked by the arrows on top of the 
plotting region) . The median was chosen because it is not 
influenced by such values. The flattening of the observa- 
tions towards the cluster center seems to be an artefact 



in the data. A comparison with newer data ( Clarke et al 
1999~p ]larke et al., 200C ) shows that the inner- most bin 
is also perfectly in agreement with the simulations. 



2.2. Initial conditions 

For an accurate calculation of the synchrotron emission 
the resolution of the magnetic field is crucial. We therefore 
perform simulations in a CDM Einstein-de-Sitter cosmol- 
ogy (n° m = 1.0, fl A = 0, H = SOkms^Mpc- 1 ), with 
twice the mass resolution of the simulation described in 
Dolag (1999). We simulate a set of ten different realisa- 
tions which result in clusters of different final masses and 
different dynamical states at redshift z = 0. 

As shown in Dolag 1999, the initial field structure is 
irrelevant because the final field structure is mostly deter- 
mined by the dynamics of the cluster collapse. The initial 
field strengths in this simulations where chosen so that the 
final field strength reproduces th e obs erved Faraday rota- 
tion data collected by Kim et al. ( 1991 ). The central values 



of the ten clusters vary between 0.5/LtG and 2/l/G, whereas 
Sim. 9 reaches values of rs 1/xG. The simulated magnetic 
field strength drops by two orders of magnitudes from the 
central part out to the virial radius (see Fig. Q). 

Given the large uncertainties a field strength twice as 
high may also be consistent with the Faraday data. There- 
fore we also give results for the extrapolated case of field 
strength twice as high. We note that neither the simula- 
tion, nor the Faraday rotation measurements can exclude 
the existence of stronger small-scale (< 1 kpc) magnetic 
field components. Here, we assume only the presence of 
large-scale (> 20 kpc) and relatively weak fields. A com- 
parison of the synthetic Faraday rotation measurements 
with the observations is shown in Fig. ||. 

3. Hadronic Secondary Electron Injection 

A CRp population may be well described by a power-law 
in kinetic energy E p 



GeV VGeVy 



dE p dV . 



(1) 



We assume a lower cutoff at E Pim - m 1 GeV. The con- 
tribution of He (assuming a 24% mass fraction) to the 
cosmic rays and thermal gas is implicitly taken into ac- 
count in our calculations by replacing proton densities by 
the correct nucleon densities without further notification. 
We choose the normalization N p (x) so that the (kinetic) 
CRp energy density ecRp( £C ) is a constant fraction of the 
thermal energy density Sth(x) of the ICM. 



£CRp = ^CRp £th 



N p (x) (E, 



^p,mm 

GeV 



2-a v 



GeV 



(2) 



Since the spectral index is sufficiently steep in our model 
(a p w 2.4) the upper cutoff was set to infinity. 

The CRp interact hadronically with the background 
gas and produce pions. The charged pions decay into sec- 
ondary electrons (and neutrinos) and the neutral pions 
into gamma-rays: 



7T° -> 2 7 



The resulting electron injection spectrum 
q e (x, E e ) dE e dV peaks at the energy of m % c 2 /Am 35 MeV 
due to rough energy equipartition between electrons and 
neutrinos in the charged pion decay (Mannheim & 
[jchlickeiscr, 1994 ). The radio emission is produced by 
electrons with an energy of several GeV (depending on 
observing frequency and magnetic field strength), which 
have to result from protons of even higher energies. For 
this energy range the electron production spectrum can 



be calculated following Mannheim & Schlickciscr (1994) 
to be 



q e (x,E e ) 



I a 



pp 



cn p {x) 



N p (x) / 24£ c 
GeV ^ GeV 



(3) 
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where a v 



32 mbarn is the inelastic p-p cross section, 



and rip is the target proton density . We note, that the neu- 
trino production spectrum is q v (x,E v ) = 2>q e {x,E v ), and 
the 7r°-decay induced gamma-ray spectrum is q 7 (x, E-y) = 
3q e (x,E 7 /i)/(16a p - 8) for E 1 »m # c 2 /2 « 68 MeV. 

The steady-state CRe spectrum is shaped by the in- 
jection and cooling process and is governed by 



„ 3 



ix. E, )Ux.E, 



q e (x,E c ) . 



For E e (x,p) < this equation is solved by 
1 



fc{x,E c 



\E e (x,E e )\ J E , 



dE' c q c (x,E' e ). 



(4) 



(5) 



„ 5 



The cooling of the radio emitting CRe is dominated by 
synchrotron and inverse Compton losses giving 



E e (x,E c ) = 



4<t t c ( B 2 (x) 



3 m 2 c 4 



d2 
CMB 

8tt 



Ei 



(G) 



where B(x) is the local magnetic field strength and 
£?Q MB /(8 7r) is the energy density of the cosmic microwave 
background expressed by an equivalent field strength. The 
CRe have therefore a power-law spectrum 



JeK ' cJ GeV VGeVy 



(7) 



with a v 



N e (x) 



3 tt P 

2 7 3 2 tt 24-t( Q >=- 
4a„ - 5 



i and 



cr T GeV B 2 (x) 



n p (x) N p (x) 
B 2 

CMB 



.(8) 



The radio emissivity j v at frequency v and per steradian 
is given by 



j v = c 2 (a c ) N c (x) Bj_(x)~ 



(Lang, 1999) with ci 



v 

2 /Co _™3 



C 2 (Q! c ) = 



3eGeV7(2 7rmf c 5 ) , 
3a r — 1 \ _ / 3a r 



167T m c c 2 a c + 1 



7 

3p 



12 



12 



(9) 



,(10) 



and B±(x) being the magnetic field component within 
the image plane. Integration of j v (x) along the z-direction 
gives the surface brightness of the radio halo. 

Thus we find that the radio emissivity as a function of 
position in the hadronic secondary electron model scales 
with the ICM properties as 



X, 



Bl» +1 
CR P n P £ th 2 2 

D CMB 



(11) 



where a„ = (a e — l)/2 = (2a p - l)/3. It follows that 
the emissivity of strong field regions (B 2 3> Bq MB ) is 
nearly independent of the magnetic field strength, whereas 
in weak field regions (B 2 <C S^, MB ) it is very sensitive to 




Fig. 3. The upper left panel shows the Effelsberg radio 
map of the radio halo in the Coma cluster from Deiss et 
al. (1997). The central radio halo, the peripheral cluster 
radio relic 1253+275, and the so called radio bridge be- 
tween them are visible. The other panels show synthetic 
radio maps of a simulated cluster (Sim. 9) at 1.4GHz 
in three spatial projections. For brtter comparison they 
are smoothed to the resolution of the radio observation. 
Contour levels are 0.1, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4 
m Jy / arcmin 2 . 

it. The factor n p e t h is nearly identical to the X-ray emis- 



sivity l x 



(kTe) 1 / 2 so that for a radially decreasing 
magnetic field strength one always gets a radio luminosity 
which declines more rapidly than the X-ray emission: 



3v 
lx 



- — <~ A-i 



CRp 



(*T e ) 



1/2. 



Bt 



(12) 



R2 I R2 ' 
D CMB 

The normalization constant of Eq. |ll| depends on the as- 
sumed scaling between CRp and thermal energy density. 
We fix this relation by comparing Coma to a similar simu- 
lated cluster (Sim. 9) and assume it to hold for all clusters 
in order to predict their radio luminosity. 

4. Synthetic Radio Maps 

To obtain the synthetic radio maps we calculate all quan- 
tities needed in a 5Mpc 3 box containing 200 3 grid cells, 
centered on the cluster position and distribute the CRe 
as described before. From this data cube we calculate the 
maps by projecting the emission along the three spatial 
directions. As the radio emission depends on the angle be- 
tween the magnetic field and the observing direction, this 
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Fig. 4. The radial profile (centered on the emission 
weighted maximum) of the Coma cluster (thick line) is 
compared to the profiles obtained from the three projec- 
tions of the simulation shown in Fig. ^ (thin lines). The 
dotted line represents the mean profile over all ten simu- 
lated clusters. The horizontal line marks a crude estimate 
of the noise enclosed in the observations. 



leads to different total radio luminosities in different ob- 
serving directions. As the magnetic field in the simulations 
is nearly isotropic, the differences are small. Additionally 
projection effects of substructure in the underlying density 
distribution can lead to different radial profiles as seen in 
Fig. | 

As a reference point we chose one simulated cluster 
(Sim. 9), which has an emission weighted temperature (9.4 
keV) and a velocity dispersion (1200 km/s derived from 
the virial mass and radius) comparable to the values ob- 
served in the Coma cl uster of galaxies (9 keV (Donelly et 
al, 1999), 1140 km/s ( [Zabludoff et al., 1993|)). This clu s- 



ter i s in a post-merger stadium , like Coma ( purns et al 
1994 |; [$ravo-Alfaro et al., 2000| ) . 

We assume a p = 2.375, which gives a v = 1.25, 
typical for cluster radio halos, and adjust A"cr p to 
0.14 (i?p jln i n /GeV)" ' 375 so that this cluster reproduces 
the observed total 1.4 GHz radio emission of Coma. In 
the case of a twice the field strength Xcr p needs to be 
0.04(.Ep im i n /GeV)~ a375 . Note that in this case the mag- 
netic fields and the CRp are in equipartition for -Ep^im = 
1 GeV. 

The synthetic radio maps are shown in Fig. [| in com- 
parison with the radio map of Coma. In the calculation 



of the radial radio emission profile of Coma the so called 
radio bridge of Coma (see Fig. [|, top-left panel at x = 3.5 
Mpc, y — 1.5 Mpc) is included. This can be justified be- 
cause a similar bridge is also visible in the X-ray emission, 
which suggests that the radio bridge results from ICM sub- 
structure. This substructure is believed to be the core of 
a smaller merging cluster which already passed the core 



region of Coma (Burns et al., 1994; Bravo- Alfaro et al 



2000). The radio relic in Coma (see Fig. |3|, top-left panel 
at x = 4.5 Mpc, y = 1 Mpc) does not affect the compari- 
son of the radio profiles due to its peripheral location. 

Although we cannot expect this simulated cluster to 
match the morphology of the radio halo of Coma exactly, 
the radial emission profiles in the y- and z-projection are 
quite comparable (Fig. ^). In the x-projection the two sub- 
lumps are projected on top of each other. This leads to a 
steeper radial radio profile (Fig. Q). The second simulated 
cluster with flat radio emission profile is also in a post- 
merger stadium with substantial substructure. 

The large extension of low level radio emission of Coma 
is not reproduced by the model. As the simulation has its 
highest resolution in the center, the simulation might un- 
derestimate the small scale magnetic fields, and therefore 
the radio emission, in the outer regions. This could explain 
the less extended synthetic radio halo. Another possibil- 
ity is that e.g. diffusion established a more extended CRp 
profile than it is assumed in this work. 

The observed radio halos do not show any radio polar- 
ization. We estimate for our simulated field configuration 
the expected degree of radio polarization following the 
treatment of internal Faraday rotation of Burn (1966). 
The synthetic polarization maps has peak values of the 
order of 10 -6 . This undetectable signal would be further 
reduced by beam depolarization due to the small-scale po- 
larization pattern. 

It is difficult to compare morphological different clus- 
ters of galaxies. But Eq. [l2] indicates a close relation be- 
tween radio and X-ray emissivities within the hadronic 
halo model. This should be approximatively conserved in 
the projected fluxes, which allows for a less morphology- 
dependent comparison of different real or simulated clus- 
ters. For Sim. 9 we compare X-ray and radio surface 
brightness point by point. In order to process the data 
in the same way as real observational data we average 
over 250 kpc boxes and use the RMS as an error estimate 
(plus systematic errors of 5Lx = 10 -7 ergs -1 cm -2 and 
<5L„ = 10 -21 ergs -1 cm -2 Hz -1 ). The results are given in 
Fig. H For instrumentally detectable regions (L v > bL v , 
Lx > SLx ) the radio and X-ray emissivity can be related 
via a power law L u = aL x , where c 



2.23 • 10 -13 and 
ive for four 



b = 1.26 in cgs units. Govoni et al. (2000) 
clusters a point by point radio to X-ray count rate correla- 
tion. Their values (b = 0.64±0.07for Coma, b = 0.82±0.04 
for A2319, b = 0.94 ± 0.04 for A2255, b = 0.99 ± 0.05 for 
A2744) are significantly smaller, which is also supported 
by the visual impression of the more concentrated simu- 
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lated halo in comparison to Coma (Fig. ||). This indicates 
that the radial CRp profile should be more extended than 
the thermal energy of the ICM, as we have assumed in 
this work, if radio halos are of hadronic origin. 

5. Gamma- Rays and Neutrinos 

The luminosity of the simulated cluster Sim. 9 in neutrinos 
and gamma-rays is Q V ^(E V ) = 5.5 • 10 45 E~ 2 b /Ge\- l b 
and Q 7 (£ 7 ) = 5.9 • 10 45 K^/GeV" 1 ' 5 . This leads to 
a gamma ray flux of < 5.2 • 10 _8 cm _2 s _1 above 100 
MeV (< since the spectrum is significant below the 
above power law near E 1 = 100 MeV). This is compa- 
rable to the upper limit 4 • 10~ 8 cm _2 s _1 for the flux 
from Coma measured by EGRET (Sreekumar et al., 
1996). In the case of twice the field strength one finds 
Q V ^{E V ) = 1.5 • lO^-^/GeV- 15 , and Q 7 (£ 7 ) - 
1.6 ■ 10 45 £- 2,5 /GeV -1 - 5 , which leads to a flux above 100 
MeV of < 1.4- 10~ 8 cm" 2 s" 1 well below the EGRET limit. 
For comparison with gamma ray and neutrino flu x pre - 
dictions for cl uster s of galaxies see Volk et al. (1996) ), 
Enfilin et a l. ( |l997|) , Colafrancesco & Blasi fll99s| ), and 
Blasi ( |1999D . 



6. The Temperature-Radio-Luminosity Relation 

Clusters with radio halos exhibit a strong correlation be- 
tween the radio luminosity (at 1.4 GHz: Li.4GHz) and 
the ICM temperature (Tx) flLiang, 199£ ; Colafrancesco, 
1999] ) . This is expected because of the known correlations 



of cluster size and energy content with temperature (Mohr 
fc Evrard, i997j ; |Schindler, 1999| ; |Mohr et al., 1999| ). But 
the shape of the Tx~E\agy^ relation is an important 
touchstone of every radio halo theory. It can be seen in 
Fig. H that the hadronic halo model reproduces the ob- 



Fig. 6. Temperature-radio luminosity relation. Th e dat a 
points with the error bars are taken from Liang ( 1999 ), 
the diamonds are the simulated clusters. Bigger symbols 
mark Coma and Sim. 9. 



served relation if a cluster independent £cR P /£th ratio is 
assumed. 



7. Conclusion 

A set of simulated clusters of galaxies in a CDM EdS 
cosmology, including fully MHD simulated magnetic field 
structures, was used to test a conceptual simple hadronic 
secondary electron model for radio halo formation of clus- 
ters. The seed magnetic fields were chosen so that the final 
magnetic fields reproduce the observed Faraday rotation 
measurements. It was assumed further, that in all clusters 
a spatially constant ratio between the cosmic ray proton 
and the thermal energy density exists. We normalized this 
ratio by comparison of a simulated radio halo of a Coma 
like cluster in our sample to the observed radio halo of the 
Coma cluster. With the further simplifying assumption of 
an unique power law proton energy spectrum, the result- 
ing radio halo, gamma ray and neutrino emission could be 
calculated. The main results are: 

1. A cosmic ray proton to thermal gas energy ratio of 
4 ... 14 % (£p,min/GeV)~ ' 375 seems to be required 
to reproduce the total flux at 1.4 GHz of the Coma 
cluster with the simulated Coma like cluster. -E^mm 
is the lower kinetic energy cutoff of the protons with 
assumed spectral index a p « 2.375. 

2. The radial profile of the simulated Coma like cluster, 
is comparable to that of Coma. 
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3. But a detailed point by point radio to X-ray compari- 
son, reveals a relation (£i.4GHz ~ L\, b = 1.26) which 
seems to be steeper than the observed relation (b < 1). 
This would imply that if radio halos are of hadronic 
origin, the parent cosmic ray proton population needs 
to have a spatially wider energy density distribution 
than the thermal gas. 

4. The expected radio polarization is of the order of 10~ 6 . 
This is mainly a property of the magnetic field con- 
figuration, and should also hold for other radio halo 
models. 

5. The expected gamma ray fluxes are below present ob- 
servational limits, but possibly in the range of future 
instruments. 

6. The set of 10 simulated clusters reproduces the ob- 
served radio luminosity-temperature relation of clus- 
ters of galaxies surprisingly well. 
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